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It is well-known that the quantum tunneling makes conventional perturbation series non-Borel summable. We 
use this fact reversely and attempt to extract the decay width of the false- vacuum from the actual perturbation 
series of the vacuum energy density (vacuum bubble diagrams). It is confirmed that, at least in quantum me- 
chanical examples, our proposal provides a complimentary approach to the the conventional instanton calculus in 
the strong coupling region. 



1. Introduction 

Usually the quantum tunneling is regarded as 
a purely non-perturbative phenomenon. This 
widely accepted picture emerges from the obser- 
vation like this: Consider a simple Hamiltonian 
of a metastable system 
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Since the potential energy is not bounded from 
below, the quantum system should be defined 
by the analytic continuation from g < 0. This 
continuation produces an imaginary part of the 
quasi-ground state energy eigenvalue, that is 
physically interpreted as the total decay width of 
the quasi-ground state due to quantum tunnel- 
ing. In the WKB (or instanton) approximation, 
we find 
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where 5*0 is the instanton 111] action and 5o = 8 
in this model. Since eq. (g) vanishes to all or- 
ders of the expansion on the tunneling effect 
is invisible in a simple expansion with respect to 
the coupling constant. In other words, a sim- 
ple (truncated) sum of the conventional Rayleigh- 
Schrodinger perturbation series of the ground 
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state energy, 
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where cq = 1/2, ci = -1/32, ca = -7/1536, 
• • •, is always real for g real. Therefore, one has 
to utilize a certain non-perturbative technique to 
evaluate the tunneling amplitude (|^): This is the 
conventional wisdom. 

However, it is interesting to note that extensive 
studies started in the seventies have revealed 
the close relationship between the tunneling am- 
plitude in the weak coupling region (|^) and the 
large order behavior of the perturbation series (||) . 
In fact, it is possible to show that P|, 
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In general, at least in super-renormalizable theo- 
ries, one finds the following correspondence: 

Tunneling rate in the weak coupling (5) 
^ Higher order perturbation coefficients. 

This relation is interesting because the right hand 
side is hard to be evaluated, while the left hand 
side is tractable. The large order behavior of per- 
turbation series in various systems has been in- 
vestigated on the basis of the connection (||) ||] . 

Furthermore it can also be shown that the 
0(17) correction in (||) gives rise to the 0(l/n) cor- 
rection for the large order behavior (^. There- 
fore, by extrapolating the relation (|5|), we expect 
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the following: 

Tunneling rate in the strong coupling (6) 
<^ Lower order perturbation coefficients. 

This possibility seems astonishing, because, al- 
though it is very difficult to systematically eval- 
uate the left hand side, we can certainly com- 
pute the right hand side: This is the basic idea 
behind our approach |3,^. (A similar approach 
was first proposed in |5| on the basis of another 
kind of resummation method.) We stress that 
the higher order corrections to the leading in- 
stanton approximation d) is difiicult: The difii- 
culty is not only technical but even a principal 
one — one has to resolve the "mixing" between 
non-trivial configurations and perturbative fiuc- 
tuations to avoid the double-counting. Also such 
an instanton expansion, i.e., an expansion with 
respect to exp(— 5*0/3), is expected to be asymp- 
totic at best. Therefore if a certain method can 
bypass these complications of the instanton cal- 
culus, it will provide a complementary approach 
to the tunneling phenomenon in the strong cou- 
pling region. This is what we want to propose in 
the next section. 



2. Borel resummation of vacuum bubbles 

As the natural extension of (^, we consider the 
0{N) symmetric cj)'^ model whose action is given 

by 
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where 0^ = 4> ■ 4>'- Note that the system is 
metastable for 5 > 0. We consider the super- 
renormalizable cases [D — 1, 2, and 3) and as- 
sume the appropriate renormalization. Then the 
standard instanton calculus gives the follow- 
ing expression for the imaginary part of the false- 
vacuum energy density. 
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where Sq denotes the instanton action, 

[-8 for £> = 1, 

So = < 35.10269 for D = 2, (9) 
I 113.38351 for £» 3, 

and An and Cd,n are some numbers arising from 
the Gaussian integration around the instanton 
and the collective coordinate integrations 

As noted in Introduction, the leading instan- 
ton approximation (H) is reliable only for g <C 1. 
Therefore, following the basic idea (^, we start 
with the conventional perturbative expansion of 
the vacuum energy density, i.e., a sum of the vac- 
uum bubble diagrams 



(10) 



ri=0 



Once the perturbation series ( p^ ) is obtained, we 
construct the Borel-Leroy transform: 
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Then we define the vacuum energy density by the 
Borel integral 
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This representation is justified by the Borel 
summability of 5 < case M combined with sev- 
eral plausible assumptions [Q . 

Now, the fact that all the perturbative coeffi- 
cients in (except cq) have the same signs is 
important (see (^ and (^). The perturbation se- 
ries of non-alternating sign is usually called non- 
Borel summable because the Borel transform ( |TT| ) 
develops singularities on the positive real axis, 
and thus the integration (12) is ill-defined with- 
out the e-separation. However, these singularities 
are welcome for us because these are precisely 
the points where the imaginary part emerges — 
Borel singularities detect the quantum tunneling. 
In fact, the leading instanton approximation (^) 
for (? <C 1 implies the nearest singularity from the 
origin is a square-root branch point 
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When substituted in (jlj), the branch cut repro- 
duces the imaginary part (||) by choosing the up- 
per contour. 

Thus, in principle, we construct the Borel 
transform ( ITTl ) from the actual perturbation coef- 
ficients c„ and substitute it in ( |l^ to extract the 
"non-perturbative" information. However, this 
way of working is impossible in practice because 
the radius of convergence of the series (|l^) is fi- 
nite (= So) due to the singularity (13). In order 
to perform the Borel integration (12), we have 
to analytically continue the series ^ll|) outside 
the convergence circle: This is impossible with- 
out knowing all the perturbative coefficients. For- 
tunately, as is well-known this difficulty can 
be avoided by the conformal mapping trick. To 
apply this trick, the knowledge on the position 
of the nearest singularity (p^, i.e., the instanton 
action 5*0, is important. 

In this way, we obtain the Pth order approxi- 
mation of the imaginary part 
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Note that eq. ( |l^ is solely expressed by the first 
P perturbative coefficients c„ and the value of 
the instanton action S'o (^. In this integration, 
the contribution around the origin = is pro- 
portional to exp(— 5*0/(7), the leading instanton 
behavior. Other parts might be regarded as the 
higher order corrections to it: This is true at least 
in quantum mechanical cases |^,^. 

3. Numerical tests 

The validity of our formula (^4|) has been tested 
numerically in quantum mechanics D = 1 and 



D ~ 2, tunneling on line ||,^. For D — \, per- 
turbative coefficients of the vacuum energy, i.e., 
the ground state energy, to very high orders are 
available. The first several coefficients are 
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With the aid of computer, it is not difficult to 
compute c„ to, say, n = 50. The exact com- 
plex quasi-ground state energy is also available 
by a numerical diagonalization of the Hamilto- 
nian. Therefore, we can compare our formula (|lj) 
and the instanton result (||) with the exact value 
of the imaginary part. For the detail, we refer 
the reader to Refs. For D — \, we see an 



excellent convergence of (14) to the exact value 
in a whole range of the coupling constant; the 
proposal in fact gives rise to the improvement of 
instanton calculus, as we announced. For exam- 
ple, for (7 = 4, which belongs to the strong cou- 
phng region eq. (^ with P = 5 (thus only 
the first five perturbative coefficients!) gives only 
a few percent error, while the leading instanton 
approximation (^ is about two times larger than 
the correct value. 

For D — 2, the situation is yet not obvious. 
We have calculated vacuum bubble diagrams to 
five loop orders under the same renormalization 
condition assumed in (^. We found (with an ap- 
propriate shift of the origin of the vacuum en- 
ergy) B: 
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Unfortunately, any converging behavior could not 
be observed and it was impossible to draw a def- 
inite conclusion on the convergence of (|lj) . It 
is not clear whether this is due to the lack of or- 
ders of the perturbation series, or there exists a 
fundamental obstruction for our approach which 
we did not encounter in quantum mechanics. To 
clarify this point, much higher order perturbative 
calculation is certainly desirable. 

4. Discussion 

We have proposed a new approach to the tun- 
neling phenomenon in super-renormalizable field 
theories. Our approach utilizes only the infor- 
mation of the conventional perturbation series 
around the naive vacuum (f> = (and the value 
of the instanton action). We have verified nu- 
merically that, at least in quantum mechanical 
cases (D = 1), we can extract a very accurate tun- 
neling rate from the perturbation series. The pro- 
cedure thus provides a complimentary approach 
to the instanton calculus in the strong coupling 
region. In I? = 2, unfortunately, the number of 
orders of the perturbation series we computed is 
yet insufficient to draw a definite conclusion on 
the convergence. 

As another test of our proposal, one of us |^ 
recently applied it to the Gaussian propagator 
model, in which the perturbative calculations in 
D = 1, 2, 3 and 4, up to the ninth order of the 
loop expansion are known It was found that 
eq. ( p^ rapidly converges in the whole range of 
the coupling constant for lower dimensional cases, 
D — I and 2, and, at least in the strong coupling 
region for higher dimensional cases, D — 3 and 4. 

As a general question, one might ask "To 
which extent we can expect the applicability of 
such a perturbative approach to tunneling phe- 
nomenon?" For example, we see that our ap- 
proach cannot work for the double-well potential 
problem: The perturbation series around one of 
potential minima is non-Borel summable and 
thus the Borel integral produces the imaginary 
part. However this imaginary part must be fake 
because the potential energy is bounded from be- 



low (i.e., no decaying process in involved). This 
fake imaginary part is believed to be canceled by 
the multi-instanton contributions ||^. In other 
words, in this system, the information of the per- 
turbation series around the naive vacuum is not 
sufficient to specify the original physical quanti- 
ties. Therefore the natural answer to the above 
question is: Such a perturbative approach to tun- 
neling phenomena is expected to be applicable 
if the system can be defined by an analytic con- 
tinuation of another system in which the Borel 
summability (of the energy density) is guaranteed. 
Eq. (^ with D = 1,2 and 3 is expected to be pre- 
cisely such a system. 

Stated differently, for our approach to be work- 
able, it is important that there is no non-trivial 
topological sector in the configuration space. The 
instanton in the system (^) is "non-topological" — 
the so-called bounce solution — and it can mix 
with the trivial perturbative sector. Our ap- 
proach systematically counts the contribution of 
the perturbative fluctuations while avoiding the 
double counting of the bounce configuration. If 
there exist a non-trivial topological sector which 
does not mix with the trivial perturbative sector, 
such as the instanton in the double well poten- 
tial, we will have to supplement the "truly" non- 
perturbative information. This is the principal 
limitation of our approach. 

When one considers a renormalizablc field the- 
ory, such as (0) with D — 4, new difficulty may 
arise The renormalon — another known 

source of the Borel singularity — emerges in gen- 
eral. The renormalon has two effects: i) The con- 
tribution of renormalons to the Borel singularity 
produces new imaginary part besides the quan- 
tum tunneling, ii) The UV renormalon suggests 
the triviality of the model. In (^ with Z? = 4, 
when g > 0, there is no renormalon singularity 
on the positive real z axis, because the model is 
asymptotically free. However the model with g < 
has the UV renormalon and is expected to be 
trivial. Therefore the meaning of the true tun- 
neling amplitude, which should be defined by the 
analytic continuation from g < 0, is not obvious. 
Presumably, in renormalizablc field theories, our 
proposal is workable only with an UV cutoff. 

Finally, we would like to comment that our 



approach is applicable not only to the quantum 
tunneling problem, but also to metastable prob- 
lems in statistical mechanics. For example, the 
free energy of the Ising model below the critical 
temperature acquires the imaginary part when 
it is analytically continued from the positive H 
to the negative H {H is the external magnetic 
field and we assumed the Ising spins are origi- 
nally aligned to the positive direction) . The imag- 
inary part might be physically interpreted as the 
inverse of the relaxation time of the metastable 
state. Our approach is applicable to the pertur- 
bation series on H, which can be computed from 
the low temperature expansion of the free energy. 
Then we expect the validity of our approach in 
the strong H which will provide a complimentary 
approach to the conventional droplet calculation. 
A study along this line is in progress. 
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